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Abstract 

We present results for masses of spin- 1/2 and spin-3/2 double-charm baryons in 
quenched lattice QCD, from an exploratory study using a non-perturbatively im- 
proved clover action at 13 = 6.2. We have studied local operators and we observe, 
after appropriate projections, a good signal for the ground states. We also present 
results for single-charmed baryons and spin-splittings for both double- and single- 
charmed states. 



1 Introduction 



Recently SELEX, the charm hadroproduction experiment at Fermilab, has reported a 
narrow state at 3519 ± 1 ± 5MeV. This state decays into A+K~7r"'", consistent with the 
weak decay of the doubly charmed baryon S+ [1] . This is the first observation of a baryon 
containing two charm quarks. BaBar and Belle at the SLAC and KEK 6-factories may also 
provide further evidence of doubly charmed baryons. 

Double charmed baryons combine the opposites of the slow relative motion of two heavy 
quarks with the fast motion of a light quark. They provide scope for testing ideas developed 
for single charm physics, such as the predicted hierarchies in lifetimes and semi-leptonic 
branching ratios and give us more room to explore predictions of exotic tetra- and penta- 
quark states (see the review by Richard [2] for further details). Unresolved issues regarding 
the recent SELEX observation itself include the fact that the observed lifetime of less than 
30 fs is much less than predicted by quark models [3] and that the observed isospin splitting 
m{ccd) — m{ccu) is rather large, about 60 times that for nucleons. The latter is, perhaps, 
the cause for publication of only one of the three states (the others being around 3460 and 
3780 MeV) observed in the first instance (see e.g. [4]). Therefore it is important to study 
these baryons further. 

The first prediction for the masses of these double- charmed baryons comes from the 
early work of De Rujula et al. [5], with later calculations from quark models and QCD sum 
rules [6-15]. 

Lattice QCD provides a method of calculating the masses of these baryons from first 
principles in a model-independent and non-pcrturbativc manner. It is interesting to com- 
pare the results from different lattice calculational techniques. Previous lattice calculations 
have used the D234 action [16] and NRQCD [17]. NRQCD is less suitable for charm quarks 
than for beauty quarks, and furthermore charm quark masses are very accessible to lattice 
simulation without using an effective theory. In this calculation we use a non-perturbatively 
0(a)-improved clover action [18]. Thus our results have discretisation errors (9(m^a^), with 
mla? about 25% in our simulation. Interestingly, although our double-charm baryons have 
masses larger than one in lattice units, the physical mass looks consistent for the one state 
for which an experimental number is available. 

We also study spin-splittings for charmed baryons and mesons, where the leading charm 
quark mass dependence cancels. Recent calculations using the 0{a) non-pcrturbativcly 
improved clover action find vector-pseudoscalar meson splittings in better agreement [19, 
20] with experiment than earher calculations using less- improved clover actions [21]. We 
confirm this here. For the single-charmed baryons, calculations with a tree-level clover 
action had difficulty reproducing the experimental splittings [22], while simulations using 
the D234 [16] or NRQCD [17] actions were compatible with experiment. We too find 
compatibility. For the doubly- charmed baryons, experimental data are not yet available. 
However, our results are compatible with those from the D234 and NRQCD actions. Since 
the hyperfine splitting is sensitive to the chromomagnetic moment term in the improved 
clover fermion action, we believe this shows the importance of using the non-perturbative 
value for its coefficient (csw)- A similar observation was made concerning the coupling 
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"c 
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CSS 2698(3) 




szi = l, J^ = 3/2+ 


K 


CUM, cud, cdd 2518(2) 




CMS, cds 2646(2) 




CSS 



Table 1. Summary of charmed baryons. Valence quark content and spin-parity are shown. 
The quantities Scc and su are the total spin of the charm and light quark pair respectively. The 
experimental values arc from ref. [23], averaged over isospin multiplets. The See mass is from the 
recent observation of the S+(cc(i) [1]. 

with the chromomagnetic field in the NRQCD action [17] (c4 in eq. (A5) in [17]). 
Our final results for the double-charm masses and splittings are 

See = 3549(13) (19) (92) MeV fi^c = 3663(11) (17) (95) MeV 

H*e = 3641(18) (8) (95) MeV fi^^ = 3734(14) (8) (97) MeV (1) 

H^e-Sec = 87(13) (13) (2) MeV Kc-^c = 67(9) (13) (2) MeV 

The splittings are determined by fitting ratios of correlators, which gives smaller errors 
compared to taking a difference of separately-fitted masses. Results for charmed meson 
splittings and single-charm baryon masses and splittings are in the body of the paper. 

In this paper the theoretical input and computational details of the simulation are given 
in Sections 2 and 3, whereas the results are presented and analysed in Section 4. There is 
a brief conclusion in Section 5. 

2 Baryon states and interpolating operators 

The double and single charmed baryons expected in QCD arc summarised in tab. 1. 

On the lattice, the masses of these hadrons can be calculated in the usual way from 
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the large time behaviour of two point correlation functions 

C(t) = 5](0|J(x,t)J(0)|0) (2) 

X 

where the J's are interpolating operators with quantum numbers to create or annihilate 
the state of interest. The choice of operators is not unique. 

For the spin-1/2 doubly- heavy baryon states, a simple operator is 

J^^eahcK (h'^l^Cl"), Shh^l (3) 



where a, b, c are colour indices, C is the charge conjugation matrix and the h and / fields 
stand for generic heavy and light quarks. 

In 5"- wave baryons with two identical quarks (heavy quarks in our case) , the two quarks 
cannot couple to spin zero and the only possibility is Shh — 1 (symmetric in both spin and 
flavour). The component s^h = as well as the operator —eahc^ {h^'^ I5C h'^) vanish. The 
coupling of the light-quark spin to Shh = 1, however, can also generate the spin 3/2 states, 
^*hh and in tab. 1. 

An interpolating operator for the spin-3/2 states can be obtained by replacing 75 with 
7^^ in eq. (3). 

Jl^^eabch; (h'^^YCl"), snh^l. (4) 

This operator also couples to spin-1/2 and projections are needed to obtain the desired 
state. The spin-1/2 masses from and are equal since there is only one spin 1/2 baryon 
in the situation where two quarks are identical. We have directly verified this property in 
our simulation. 

Another operator, used for spin-3/2 double heavy baryons [17,24] is 

J!;^eabcl',{h'^rCh''). (5) 

We have also tried this operator and we see no reason to prefer one over the other. Indeed, 
both give a good overlap for the ground state and the masses extracted turn out to be 
equal as expected. 

For the operators and (or J^) the 2-point functions in eq. (2) have the following 
large-time behaviour 

^^^77 - 5^(0|J^(x,t)J^(0)|0) (6) 

^ ^1/2 (P+),, e-v.*+ Zf/, (P_)^- e-</^*. 
C2{ty^, = $^(0|j;(x,t)J^(0)|0) (7) 

X 
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Figure 1. Comparison of the effective mass plots for the double heavy operators J^, Ji^ and 
, with Kh = 0.1222 and ni = 0.1351. In each plot the upper points are for spin-3/2 and the 
lower points for spin-1/2. The spin-1/2 plateaus are the the same for Jy and (left) while both 
plateaus coincide for the Jt^ and operators (right). 



where the projection operators are defined by 



-^3/2 



1+70 
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37 7' , 



M/2 



1-70 



37 7' 



(8) 



Contributions of negative parity states are removed by projection with P+. The negative 
parity states can, in principle, be detected by using the projector P_, but in our simulation 
they are much noisier. We show an example of the signals from the operators J^, and 

in fig. 1. As stressed above, spin-1/2 masses extracted using the three operators are 

equal, while the choice between and makes no difference for the spin-3/2 mass. 
For a baryon containing a single heavy quark, a common choice of operators is 



Sll 
Sll 



0, 



1, 



(9) 
(10) 



for the states su = and su = 1 in tab. 1, respectively. In our simulation, the light quarks 
li, I2 carry different flavours but the same masses. 
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amp 


arriv 


0.1344 
0.1346 
0.1351 
0.1353 


0.300(2) 
0.276(2) 
0.210(3) 
0.177(2) 


0.397(4) 
0.383(5) 
0.352(11) 
0.340(15) 



Table 2. Light pseudoscalar and vector meson masses. The fit interval is [12 — 28]. Our time 
counting starts from 0. 

It should be noted that for baryons three different quarks, i.e., hlil2 (or //ii/i2), these two 
operators correspond to different physical spin- 1/2 states with Su = and 1 respectively, 
the latter one often being denoted by a prime. This is evident from the experimental 
masses of Sc and 5^ in tab. 1. 

3 Details of the Simulation 

Our simulation was made using the code FermiQCD [25] on a PC cluster. In this study 
100 quenched gauge configurations were generated at /5 = 6.2 on a volume of 24^ x 64 with 
1000 heatbath steps for the thermalisation followed by 200 heatbath steps to separate each 
gauge configuration. These numbers were decided upon after an autocorrelation study on 
the average plaquette values. 

Four light quark propagators around the strange quark mass and three heavy quark 
propagators around the charm were calculated using the following values of the hopping 
parameters: 

• Ki^ 0.1344 , 0.1346 , 0.1351 , 0.1353; 

• 0.1240, 0.1233, 0.1222. 

The propagators were generated by the Bistabilised Conjugate Gradient method [26] for 
the non-perturbatively improved clover action [18]. 

Since the signal is satisfactory with local interpolating operators, no smearing was 
required. The statistical errors were estimated by a jacknife procedure, removing 10 con- 
figurations at a time from the ensemble. 

3.1 Lattice spacing and quark masses 

To fix the lattice spacing, we used the method of lattice planes [27]. In other words, we 
perform the following fit to the hght vector and pseudoscalar masses in table 2, 

amy — C + L [ampY . (11) 
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Figure 2. Light vector masses as a function of squared light pseudoscalar masses. The interpo- 
lated kaon and extrapolated pion masses are also shown. 

This is shown in fig. 2. From the physical values of mx* and rriK, the inverse lattice spacing 
is found to be 

a-^ = 2.6(1) GeV. (12) 

Terms of 0{{amp)^) in eq. (11) turn out to be irrelevant and do not affect the above 
estimate (compare the linear and quadratic fits in fig. 2). For illustration, the values of the 
pseudoscalar masses in tab. 2 converted to physical units are 

mp = {779 , 716 , 546 , 459} MeV. (13) 

These span the kaon mass while the pion is instead quite far away. For this reason, we 
interpolate for the strange quark and extrapolate for the up/down masses. This is also the 
reason for using K, K* to fix the lattice spacing. 

For the heavy sector, the D^-meson mass is within our range of simulation. This is 
evident once the heavy-light pseudoscalar masses in tab. 3 are interpolated to the strange 
mass (through the lattice plane method) and expressed in physical units 

rrihs = {1.83 , 1.89 , 1.98} GeV. (14) 



7 







arrihii' 






amp 


jp 1 + 


jP 1 + 


tP 3 + 
-2 


jP 1 + 

-2 


tP 3 + 
•J -2 






Sn' — 


S^i' — 1 


sw — 1 


Shh — 1 


Shh — 1 


n 1 o/in n 1 Q/i /I 
U. iZ4U-U. io44 


n 71 c (n\ 
U. ( io(^Zj 




n no o /'^ \ 
U.yoot^O ) 


l.UUo^oj 


1 oo^ / o\ 
i.62b[o) 


1 O /I \ 

l.oo4(o j 


U.izoo-U.io44 


U. /4U (^Zj 


0.975(5) 


1.010(6) 


1.029(6) 


1.368(3) 


1.395(3) 


n 1 ooo n 1 Q /I /I 
U.izzz-U.io44 


U. / 10\Z ) 


i.UUol^D ) 


i.U44(^0 J 


i.uozi^o ) 


i.4oo(^o j 


i.4oy(^o j 


^^(^18^1-^.1044 




i.UUol^zo ) 




i.Uo/(^ozj 


i.Uool^oi j 


i.44z(^0 / ) 


U.iz4U-U.io4D 


n 71 n^'o^ 
U. / iU(^z J 


n no y( /'i^ \ 
0.yo4(^D ) 


n n7o/'7\ 

ly ( ) 


n ^^o/'7^ 

o.yyz(^/j 


1 Q 1 O /' y( \ 

i.oio(^4j 


1 O /I 7/" \ 

i.o4/(4j 


n 1 0QQ n 1 Q/1 R 


n 7QQ/'0^ 

U. / ooyZ ) 


0.9oD(Dj 


0.994(7) 


1.013(7) 


1.360(4) 


1 DO O /0\ 

1.388(3) 


U.izzz-U.io4D 


\J.I0I[Z) 




1 nOQ^7^ 

i.Uz5( / J 


1 n/lR^7^ 
i.U4D( / ) 


1 /loc;^Q^ 
i.4zo(^oj 


1 /I c:o^Q^ 
i.4oz(o j 


If n 1 Q/IR 
'«charm-U.io4D 




U.yc54(^zc5 J 


1 noQ^'Q/l^ 
i.Uzo(^o4 J 


i.U4i(^oo J 


1 /l1^J^'c;7^ 
i.4iD(^0 / J 


i.44z(^o / J 


U.iZ4U-U.iOOi 


U.Dyi l^o ) 


n Q7Q f 1 


n non^i Q^ 
\j.\)l\)y\.6) 


u.y4o(iUj 


i.Z\)iy4:) 


1 Qon^c:^ 

i.ozy^oj 


n 1 0QQ n 1 Qcii 

U.iZOO-U.iOOi 


U. / i-^yo ) 


0.900(10) 


0.951(13) 


0.966(10) 


1.339(5) 


1 o7rt / \ 

1.370(5) 


n 1 ooo n 1 Qc;i 

U.iZZZ-U.iOOi 


n 7/1 Q ^Q^ 
U. / ^oyo ) 


0.934(10) 


0.984(14) 


0.998(10) 


1.404(5) 


1.434(5) 


'^charm-0. 1351 




0.928(27) 


0.979(39) 


0.993(34) 


1.395(58) 


1.425(58) 


0.1240-0.1353 


0.683(3) 


0.854(13) 


0.903(17) 


0.915(12) 


1.287(5) 


1.322(6) 


0.1233-0.1353 


0.706(3) 


0.876(13) 


0.923(17) 


0.935(12) 


1.330(5) 


1.363(6) 


0.1222-0.1353 


0.740(3) 


0.910(13) 


0.956(18) 


0.967(12) 


1.395(6) 


1.427(6) 


'^charm-O.lSSS 




0.904(28) 


0.951(43) 


0.962(34) 


1.385(58) 


1.418(58) 



Table 3. Double and single-heavy baryon masses in lattice units, together with pseudoscalar 
masses. The fit intervals are [16 — 28] for double and [15 — 25] for single-heavy baryons. 

3.2 Analysis of the baryon masses 

Since Kcharm IS rather close to our third = 0.1222^, as the first step in our analysis, 
we interpolate the quantities of interest, viz. the single and double heavy baryon masses, 
to the charm mass. In practice, this procedure is implemented by doing for each ki the 
following fits 

arrihhi ^ Q + Li anihs , aruhii = C; L[ anihs- (15) 

Quantities at the charm mass, uicd and rricii are obtained by putting nihs — fno^- This 
interpolation is shown for the double heavy case in fig. 3. With the charm mass fixed, the 

naive linear fit in I/kh to the masses in eq. (14) gives Kcharm = 0.1224(9). 
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Figure 3. Spin-1/2 double-heavy baryon masses for all k combinations. For each we fit the 
heavy quark mass dependence using the heavy-strange pseudoscalar meson mass. The fit function 
is given in equation (15). The vertical dashed line indicates the Ds meson mass (in lattice units) 
used to fix the masses of the ccl spin-1/2 baryons. 



light quark mass dependence is studied using 

arricii = A + B {ampY , auicd = A' + B' {ampY. (16) 

This fit is shown for the spin-1/2 double charm case in fig. 4. The masses of charmed 
baryons containing strange and/or up/down quarks are obtained by the following substi- 
tutions for mp in the above equations: 

• mp = m^ for rricud, rriccu] 

• mp = rriK for m^su] 

% mp = m^^^ for mess, rriccs- 

where m^^^ = 2m\ — m^. In the second case, we suppose that SU{?)) breaking terms 
are negligible and obtain our estimate from states containing two mass-degenerate light 
quarks [24,27]. 
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Figure 4. Spin- 1/2 double charm state masses as a function of the square of the light pseu- 
doscalar masses. The values at strange and up/down masses are shown. 

4 Results and Discussion 

Here we recall our final values for the double-charm baryon masses. 

= 3549(13) (19) (92) MeV fi^c = 3663(11) (17) (95) MeV , . 

S*, = 3641(18)(8)(95) MeV n*^ = 3734(14) (8) (97) MeV ^ ^ 

The first error is statistical. The second error is systematic, estimated by combining in 
quadrature the effects of the following variations in our analysis: 

• changing the time fit-ranges — this contributes up to 35% of the quoted error; 

• using single or double exponential fits — we saw no change in our lowest state masses; 

• linear versus quadratic chiral extrapolations — in the worst case this gives three- 
quarters of the quoted error; 

• interchanging the order of light quark extrapolations and charm quark interpolation 
— this produces no change in our results. 

Only one volume and lattice spacing was studied; investigation of discretization errors, the 
continuum limit and finite volume effects are not addressed. To account for these (and 
the effects of quenching), the third quoted error is found by rescaling our masses using the 
experimental Ac mass. 
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This work [MeV] 


Expt [MeV] 


Ac 


2227(50) (57) (58) 


2285(1) 


' — 'C 


2374(34) (23) (61) 


2469(1) 


Sc 


2377(38) (84) (62) 


2452(1) 


' — 'c 


2502(26) (40) (65) 


2575(3) 




2627(16)(48)(68) 


2698(3) 




2396(42)(122)(62) 


2518(2) 


' — 'c 


2532(31)(62)(66) 


2646(2) 


n* 


2669(21) (26) (70) 





Table 4. Our estimates for the single charm baryon masses compared to experimental values. 

The Sec mass is in good agreement with the experimental value [1] 

(Scc)expt = 3519 ± 1 ± 5 MeV (18) 

Other masses are consistent with the lattice estimates using NRQCD [17] or D234 [16] 
actions. For recent estimates in quark models or QCD Sum Rules we refer the reader 
to [15] and [13] respectively. For completeness, our estimates for the single charm baryon 
masses are given in tab. 4 along with the experimental results. Values turn out to be 
compatible with previous lattice calculations [17,22,24]. It may be noted that in ref. [22], 
a perturbative value for the coefficient csw was used in the clover action. 

We now turn to the baryon and meson spin-splittings. Our results for these are given 
in tab. 5. The values are obtained either from the difference in individually fitted masses 
(labelled "Diff" in the table), or by directly fitting a ratio of correlators (labelled "Ratio" 
in the table). When using the ratio the noise starts to dominate earlier so we restrict our 
fit to a shorter time-slice window. For the baryons we find a better signal using the ratio 
method and the difference between the two approaches becomes more apparent as we move 
away from our region of simulation to lighter quarks. We use the numbers from the ratio 
as our best estimates. 

For the double-charm baryons we observe a good signal for non-zero splittings. For 
the single-charm baryons, where experimental data is available, our results are compatible. 
This distinguishes our results from earlier ones using a less-improved clover action [22] . Our 
values are also compatible with those found using the D234 [16] or NRQCD [17] actions. 
For the mesons too our results are compatible with experiment: this improved agreement 
is also found in other recent non-pcrturbatively improved clover simulations [19,20]. 

The predictions arc more precise for {Q*,Qc) and the double charm spin doublets 
(0*^,f2cc) and (S*^, Sec), where less extrapolation is needed, but experimental numbers 
are still awaited. 
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'— 'cc 


^cc 


89(15) 


87(13)(13)(2) 
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67(9)(13)(2) 








18(51) 


49(39)(12)(1) 
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' — 'c 


' — 'c 


30(33) 
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71(3) 




fie 


43(17) 


44(16)(15)(1) 




D* - 






127(14)(1)(3) 


142(2) 




Ds 




123(11)(1)(3) 


138(2) 



Table 5. Our results for the single- and double-charm mass splittings. 

5 Conclusion 

We have presented exploratory quenched lattice results for double charm baryon masses. 
These have drawn attention after the experimental observation of the first double charm 
state last year. The calculation is done with non-perturbatively (9(a)-improvcd Wilson 
fermions at /3 = 6.2 and on a large lattice. Good signals for the ground states are observed 
without recourse to smearing. In addition, we have reported the masses of single charm 
baryons. The calculated masses look quite reasonable. We see a definite signal for non-zero 
baryon and meson spin splittings. To improve our lattice calculations, a finer lattice spacing 
is necessary together with an examination of chiral logarithms in the light extrapolations 
and a simulation with dynamical quarks. 

Experimental observations of the remaining double charmed baryon states and, in par- 
ticular their spin-splittings, would allow the lattice predictions to be checked. 
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